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Abstract. We present results of a detailed chemical analysis performed on 23 main-sequence tumoff stars having 
-3.4 <[Fe/H]< -2.2, a sample selected to be highly homogeneous in T^ff and log(^). We investigate the efficiency of mix- 
ing in the early Galaxy by means of the [Mg/Fe] ratio, and find that all values lie within a total range of 0.2 dex, with a standard 
deviation about the mean of 0.06 dex, consistent with measurement errors. This implies there is little or no intrinsic scatter in 
the early ISM, as suggested also by the most recent results from high-quality VLT observations. These results are in contrast 
with inhomogeneous Galactic chemical evolution (iGCE) models adopting present supernova (SN) II yields, which predict a 
peak-to-peak scatter in [Mg/Fe] as high as 1 dex at very low metallicity, with a corresponding standard deviation of about 0.4 
dex. We propose that cooling and mixing timescales should be investigated in iGCE models to account for the apparent dis- 
agreement with present observations. The contrast between the constancy and small dispersion of [Mg/Fe] reported here and the 
quite different behaviour of [Ba/Fe] indicates, according to this interpretation, that Mg and Ba are predominantly synthesised 
in different progenitor mass ranges. 
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1. Introduction 

The chemical evolution of the Milky Way arises from the con- 
tinuous exchange of material between stars and the interstel- 
lar medium (ISM). The surfaces of long-lived, low-mass stars 
retain the composition of the ISM from their formation, so 
the early phases of the ISM are recorded in metal-poor stars. 
Metal-poor stars exhibit enhancements of oxygen and other a- 
elements relative to iron, which are explained by the differing 
yields of supernovae (SN) of type II and type la. The O comes 
from SN II (Thielemann et at. J990,0, whereas Fe is produced 
by both SN la and SN II (Nomoto et at. ,1984.) ). It is usually 
inferred that the early Galaxy must have been enriched only 
by SN II; only more recently have SN la started to pollute 
the ISM, reducing the a-enhancement at higher metallicities 
(Tinsley 1980). This implies that the chemical composition of 
metal-poor halo stars is determined exclusively by SN II events 
and subsequent ISM mixing in the early Galaxy. 

Trends of the abundances of different elements with metal- 
licity are not the only tracers of the enrichment of the Galaxy. 
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The scatter in the relative abundances of low-metallicity stars 
indicate the level of mixing of the early ISM. Whether all of the 
observed scatter is an intrinsic property of the ISM, and thus 
due to incomplete mixing in the early Galaxy, or is possibly 
due to eiTors that exceed the formal estimates of uncertainties, 
is still unclear. Nucleosynthesis models invoke various sources 
to explain the observed abundances of heavy elements (e.g. 
Travaglio et a/. l2003l l. rising the level of uncertainty in their nu- 
cleosynthesis production. However, for some lighter elements, 
such as the a-elements O and Mg, their expected modes of syn- 
thesis do not predict differences in the ejected yields for stars of 
similar mass. Any scatter observed in the abundances of these 
elements should thus reflect more directly the effects of enrich- 
ment and mixing events in the early Galaxy. 

SN II progenitors evolve on timescales of a few to tens 
of Myrs. In the earliest phase of the Galaxy there may have 
been a period of star formation during which the zones en- 
riched by SN II have not yet mixed prior to the formation of 
the next stars. In this case, the newly forming stars will have 
different chemical compositions. Audouze & Silk J 19951 1 pro- 
posed that, at very low metallicity, clouds in the ISM could 
have been polluted by a maximum of three SN II, consistent 
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with the suggestion by Ryan et al. J1991b ) that single SN were 
sufficient for the enrichment of the gas that formed stars at 

[Fe/H] 3.5. If the early ISM was dominated by local inho- 

mogeneities, as predicted by inhomogeneous Galactic chemi- 
cal evolution (iGCE) models (e.g. Ishimaru & Wanaio ri999l 
Tsujimoto & Shigevama 119991 Argast et al. 2000, Travaglio 
et al. 120011 . and the [O/Fe] and [Mg/Fe] ratios depend on SN- 
progenitor mass or metallicity, then [O/Fe] and [Mg/Fe] should 
show a scatter in the abundances trend versus metallicity, the 
surface abundances of these stars reflecting the local enrich- 
ment of the ISM. Argast finds, for example, that although the 
ISM is well mixed at [Fe/H] > -2.0, there is essentially no 
mixing at [Fe/H] < -3.0. This gives rise to progressively more 
star-to-star scatter as [Fe/H] decreases below -2.0. 

Although there have been a large number of observational 
studies, the usually required agglomeration of results from 
different authors introduces the possibility that the reported 
scatter of the a-element abundance may be due to a lack 
of internal consistency, as discussed by Norris et al. (I2001t . 
Underscoring this possibility, Magain ( I1987> found that an ho- 
mogeneous re-analysis of literature data comprising 21 stars 
over the range -3.6 <[Fe/H]< -0.9 gave a mean Mg abun- 
dance [Mg/Fe] =0.45 and a standard deviation, as a measure 
of the scatter, of 0.1 dex (scatters of the order 0.15 dex were 
reported for [Ca/Fe] and [Si/Fe]). Ten of the stars studied by 
Magain fell in the interval -3.4 <[Fe/H]< -2.2 which we study 
in the present paper A similar result was found by Nissen et al. 
<1994> for higher metallicities: [Mg/Fe]=0.41 and standard de- 
viation 0.07 dex, over the range -2.0 <[Fe/H]< - 1 .0. Recently, 
Carretta et al. ( I2002t found that the scatter in their estimate of 
intermediate mass elemental abundances was compatible with 
their estimated observational uncertainties. The most recent 
study of Cayrel et al. (2D03"), based on extremely high-quality 
VLT/UVES observations of giant stars from the HK survey, ex- 
hibit a low scatter in the abundances of most elements, down 
to metallicities as low as [Fe/H] =-4.2. These authors report a 
standard deviation about their mean [Mg/Fe] value of 0. 13 dex, 
and as low as 0.05 dex for [Cr/Fe]. This indeed supports the 
hypothesis that there is no large scatter in the early ISM. 

Other studies, e.g., Fuhrmann et al. J1995t and Mashonkina 
et al. J2003> have reported higher scatters in the alpha-element 
ratios. Even though these studies suggest the existence of stars 
with a different enrichment history (see e.g., Shiegeyama & 
Tsujimoto ( I2003I I for a possible origin of metal-poor stars with 
low [a/Fe]), their average standard deviation is still of the 
order ~ 0.1 dex, far below what is expected in the case of 
an early inhomogeneous ISM. However, among the eight ex- 
tremely metal-poor stars ([Fe/H]< -3.5) considered in the re- 
view by Norris {2003 ), two objects show [Mg/Fe] > 1.0. This 
indicates the existence of some real differences, but only for 
metallicity much lower than previously thought. 

Clearly, studies of the mixing efficiency in the early Galaxy 
are fundamental for testing inhomogeneous GCE models. This 
can be accomplished by a detailed analysis of [Mg/Fe] in halo 
stars, and is the main purpose of this work. We concentrate 
on [Mg/Fe] because of the suitability of Mg measurements 
in stars at very low metallicity, and to avoid the still debated 
uncertainties on measurements of O. Studies of the scatter 



vs. metallicity for different stellar populations could also al- 
low one to constrain the yields from SN II. Attempts to do 
this have so far been unsuccessful for individual stars (Chieffi 
& Limongi 2002 1 because of the difficulty of computing SN 
yields that match the composition of metal-poor stars. Fran9ois 
et al. ( 2004 1 analyse how SN II yields can be constrained by the 
observations on the basis of their homogeneous GCE model, 
suggesting that major revisions of SN II yields are required to 
match the elemental abundances observed at very low metallic- 
ities. 

Iron yields are particularly difficult to predict, because 
Fe production depends on the uncertainties in the mass cuts 
adopted for SN II models as well as the degree of mixing prior 
to fallback (Nomoto et al. 1 19971 1. whereas Mg is produced in 
hydrostatic burning much farther out in the star, and hence is 
subject to fewer uncertainties. As a result, Mg could also be a 
better chemical chronometer compared to Fe; its study would 
aid our understanding of the age-metallicity relation, a key con- 
straint on studies of chemical evolution and nucleosynthesis. 

The primary motivation for this study is to produce an 
independent estimate of the level of elemental scatter in the 
early Galactic ISM, as unaffected as possible by uncertainties. 
Unlike the most recent study of Cayrel et al. (2003), whose S/N 
ratio is higher ( > 200/1, compared to our 100/1-150/1), and 
who thus obtain better precision on a single measurement, here 
we focus on the homogeneity of the sample, and demonstrate 
how our homogeneous chemical analysis leads to a better con- 
straint on the observed scatter of our derived abundances. The 
aim is to avoid the introduction of star-to-star differences that 
may lead to increased and misleading eiTors. We thus concen- 
trate on halo stars that are in the same evolutionary stage, have 
only small differences in stellar parameters, and make use of 
the same Mg absorption lines for all objects. It is of fundamen- 
tal importance to understand whether the scatter we observe 
is intrinsic to the sample or due to deficiencies in the analy- 
sis. We present the results of the chemical analysis of a sample 
of main-sequence turnoff halo stars, focusing on [Fe/H] and 
[Mg/Fe]. We describe the definition of the sample and the ob- 
servational data in ^ Particular importance has been given to 
the use of techniques that produce the highest homogeneity of 
analysis. The methods adopted and the chemical analysis are 
presented in ^ Sectionl^reports our results and the determi- 
nation of their uncertainties. In Isjwe analyse the results and 
focus on their possible implications. 

2. Definition of thie sample, observations and 
photometry 

2.1. Data 

The stellar sample was originally selected by Ryan, Norris, & 
Beers C1999 - hereafter RNB) from the surveys of Schuster & 
NissenjpHl; Ryan ([T989t : Beers et al. ( fT9^ and Carney et 
al. (1994), for the purpose of studying Li in metal-poor stars. 
This sample consists of 23 main-sequence turnoff stars with a 
narrow range of stellar parameters, the target being effective 
temperatures 6100 + 50 K < Teff < 6300 + 50 K and metallic- 
ities in the range -3.5 < [Fe/H] < -2.5. The choice of turnoff 
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stars excludes the presence of subgiants, and reduces the range 
of logarithmic surface gravity to values within a few tenths of a 
dex from an expected theoretical log(^) ~ 4.0 dex. Systematic 
errors introduced in the analysis will thus affect the mean val- 
ues of our estimates, rather than introducing large star-to-star 
differences, allowing for high internal consistency in our de- 
rived results. 

Observations were made with the 3.9m Anglo-Australian 
Telescope (AAT), with the University College London echelle 
spectrograph (UCLES) - see RNB for details. The stellar spec- 
tra are of high resolving power (R ~ 40000), high signal-to- 
noise ratio (S/N), exceeding 100/1 in most cases, and (incom- 
plete) wavelength coverage over the range A x 4900-8200 A. 
Stars were observed in different epochs in 1996, 1997 and twice 
in 1998 (hereafter 1998a and 1998b). For 12 of the 23 targets, 
spectra were co-added from different epochs. At each epoch, a 
series of multiple observations was made to increase the final 
S/N obtained. 

Photometric measurements are available from a few con- 
sistent sources, for Johnson-Cousins UBVRI colours (Ryan 
I1989> . and for Stromgren uvby (Schuster & Nissen J 19881 
11989 ') and Schuster et al. (1993 , "1996 1). For uncertainties and 
interstellar reddening we adopt the estimates of RNB (and ref- 
erences therein) and calculate the reddening for the ci index as 
E(ci) = 0.20E(b - y) (Golay 1974i. As discussed below, we 
use photometric estimates extensively for the stellar parame- 
ters due to the very low number of spectral lines in the lowest- 
metallicity stars. The strengths of the Mg b lines at 5172 and 
5183 A, even at extremely low metallicity, permits Mg abun- 
dance determinations in all of our stars. 

2.2. Galactic Populations 

The very low metallicity of most of our program stars is 
prima facie evidence for their being part of the halo popu- 
lation. However, some recent studies have claimed that the 
thick disk could extend to include stars of quite low metal- 
licities, [Fe/H] ~ -2.0, though with very low frequency. See, 
for example, Allen et al. (1991j . Chiba & Beers (2000 1, one 
star of Ibukiyama & Arimoto (I2002> . and Beers et al. (20021. 
Therefore, we further investigate the likely halo membership of 
our sample stars using the criteria of Feltzing et al. (2003 1, i.e., 
adopting a 200 km s"' radius on a Toomre (UW=(\J^+Wy^^ 
vs. V) diagram as the border between the thick disk and halo. 
This limit is suggested by their data, as well as the data of 
Fuhrmann J1998> . This is a somewhat more detailed criteria of 
interpreting the actual three-dimensional motion of the stars, as 
compared to criteria where only the V component of motion is 
considered in the identification of the population of the stars 
(see e.g., Ibukiyama & Arimoto "2002 for a description of this 
second method). U,V and W components were computed by 
Ryan & Norris <1991b ) and by Carney et al. J1994> for the stars 
in their sample. Heliocentric velocity components from Ryan & 
NoiTis are transformed to the LSR frame by assuming the solar 
motion relative to the LSR is (-9, 12, 7) km s"' in the directions 
corresponding to (U, V, W) (Mihalas & Binney'1981't. The anal- 
ysis in the Toomre diagram shows that all but two of the stars in 



our sample lie outside the 0-200 km s"' zone, i.e., they are al- 
most certainly halo stars. Of these two exceptions, BDh-3°740 
has a metallicity of [Fe/H] = -2.78, surely a halo metallic- 
ity, while BD-h26°3578 has [Fe/H] = -2.22, where both es- 
timates are based on high-S/N lA resolution spectra (RNB). 
We found that Carney's estimates of UW are lower by 50 to 
150 km s ' than those of Ryan and Norris J1991b ) for the 
six stars in common. Fuhrmann et al. ( 1998|l and Ryan et al. 
(2001) found Carney et al. photometric distances to underes- 
timate the Hipparcos distances in each case for nine stars in 
common. This possibly explains the difference in the velocity 
components. A coiTection was applied to Carney's distances of 
these two stars, which moved the velocity of BD-i-3°740 into 
the halo zone, but still left BD-i-26°3578 in the thick-disk zone. 
However, our estimate of its metallicity, [Fe/H]= -2.49, makes 
it a rather unlikely thick-disk suspect, hence we conclude that 
we are analysing stars that are sampling exclusively the halo 
stellar population. 

3. Chemical Analysis 

In this section we describe the procedure followed to derive 
the abundances of chemical elements from the observed spec- 
tra, using Kurucz's WIDTH6 code (Kurucz & Furenlid 197^. 
The code assumes local thermodynamic equilibrium (LTE) for 
the formation of the spectral lines. We measure the equivalent 
width (EW) for each line we consider to be reliable on the basis 
of a selected list of lines that is described below. A model atmo- 
sphere with the appropriate stellar parameters, measured EWs, 
and the atomic data for the spectral lines are introduced as in- 
puts to the code. The code then calculates the abundances for 
each individual line by requiring the calculated EW to match 
the observed one. The final abundance of a chemical element 
is the average of the line abundances. We discuss the details of 
this analysis in the following sub-sections. 

3.1. Line selection and atomic data 

We have given priority to establishing the reliability of each 
individual spectral line, rather than to achieving a large total 
number of lines. Our initial selection is based on the lines ob- 
served in the solar spectrum (Moore et al. 1966). We reject 
lines with possible blending. We do not include lines with gf- 
values having reported uncertainties larger than 25%, and dis- 
card lines for which, for a large part of the sample, we derive 
abundances systematically different from the mean abundance 
of all the lines of a specific element. The latter can be due to 
cases in which the line sits at the very edge of a spectral order, 
or to poorly known gf data. 

Unfortunately, no single source gives laboratory gf val- 
ues for all lines, so possible systematic differences between 
atomic data sources could introduce systematic differences in 
the abundances derived from different lines. For all lines we 
adopt recent experimental or solar ^/-values and avoid theo- 
retical ones, because of the large discrepancies that have been 
found in some cases. Whenever possible, ^/-values are aver- 
aged from several reliable sources. Because of our scientific 
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aims and the paucity of lines of other species, special attention 
was addressed to Fe I, Fe II and Mg I lines. 

Fe I - For Fe I lines we used mainly ^/-values from 
O' Brian et al. J1981> and from the Oxford group (see Table 
n]for references). However, in attempting to use Fe lines with 
^/-values listed only by Fuhr et al. ( 1988 1, we found a system- 
atic difference in the derived abundances of ~ -0.1 dex. We 
thus do not include lines from Fuhr et al. alone. 

Fe II- The importance of Fe II lines is to constrain the sur- 
face gravity of the star. All but the five most metal-poor stars 
exhibit Fe II lines in our spectra. We consider Fe II lines at 
4923, 5197 and 5276 A, and use Ki-oll & Kock ( 19871 gf val- 
ues. Unfortunately, only for half of the sample do we obtain 
good measurements; many Fe II lines are distorted by noise 
or do not pass the EW threshold described below. Moreover, 
the Fe II line at 4923 A was identified only in the observations 
from 1998a, which have a larger wavelength coverage. A com- 
parison of the abundances derived from this line and from the 
other two lines do not show any particular difference. Thus, we 
believe we do not introduce any systematic error in the Fe II 
abundance for the subset of stars observed in 1998a. 

Mg I- The Mg lb triplet around 5 172 A is the most obvious 
feature in all the stellar spectra of our sample. This allows us to 
measure Mg lines even in the most metal-poor stars, where typ- 
ically only four or five lines are usable in the observed spectral 
region. We thus rely on the triplet to obtain high star-to-star 
consistency in the derived Mg abundances. We also compare 
these with abundances derived from the weak lines at 5528 and 
57 1 1 A whenever we can measure them. For Mg I we adopt 
the solar ^/-values of Fuhrmann et al. fl995 1. These are in ex- 
cellent agreement with the values of Wiese and Martin (1980 1, 
and the theoretical calculations of Chang |1990|f and Chang & 
Tang n^90) for the Mg lb lines at 5172.698 and 5183.619 A 
and for the Mg I Une 5711.092 A. For the line at 5528.42 A, 
we notice that the gf -value of -0.5 1 used by Fuhrmann et al. 
is lower than the value -0.341 of Wiese & Martin, but its adop- 
tion leads to higher internal consistency of the Mg abundances 
derived from different Unes. 

The Mg lb triplet lines are very strong, with extremely 
broad wings at solar metallicities; even in our low metallicity 
range the Une wings are still moderately strong . Fuhrmann et 
al. J1995l t analyse Mg I lines at 457 1, 4703, 4730 and the triplet 
at 5172 mA for 56 metal-poor stars. They report no large sys- 
tematic differences in the Mg abundances derived from these 
strong lines and from weak Mg lines, nor from lines with dif- 
ferent excitation potentials. However, for the three stars where 
they identify our subset of weak lines (5528 and 571 1 A ), they 
find abundances derived from the latter to be about 0.05 dex 
smaller. In our results there is no systematic difference affecting 
all the stars. Single cases in which we find a difference between 
weak and strong line Mg abundances are more likely related to 
the uncertainties in the model parameters (as discussed below). 

The use of the two strong Mg lb lines requires a more care- 
ful investigation of the damping constant. The classical treat- 
ment of the van der Waal's broadening factor by Unsold (1955 ) 
has recently been revised. The formalism of Anstee & O'Mara 
([T998'l leads to 7ao = 2.32 X yunsoid, while the value derived 
from a fit to the Mg lb lines in the spectrum of the Sun is 



70 = 2.04 X y Unsold- This suggests a ~ 10% uncertainty in the 
damping factor, which corresponds to an uncertainty of the or- 
der of 0.02 dex in our Mg abundances, since our Mg lines lie 
at the beginning of the flat part of the curve of growth. We as- 
sume an intermediate damping constant of 2.2 x yunsoU- We 
note that the suggestion of Gratton & Sneden ( I1994> of adopt- 
ing a factor SxyunsoU would yield a worse agreement of the Mg 
abundances derived from strong and weak lines in our sample. 
We do not investigate this possibility any further in this study, 
and rely instead on the estimates above. 

This selection procedure leads to a total of 52 reliable lines 
(including 4 Mg, 32 Fe I and 3 Fe II lines), as shown in Tabled 
Only eight of these lines (of which three are Mg lines) are mea- 
surable in our most metal-deficient stars. 

3.2. Equivalent widths 

We used IRAF to continuum fit the individual orders of the 
echelle spectra, and to shift the spectra to the rest frame. 
Avoiding lines distorted by excessive noise or by cosmic rays, 
the equivalent width (EW) of each line from our list was mea- 
sured by fitting a Gaussian profile (for lines with EWs up to 
~ lOOA ) and a Voigt profile for stronger lines (i.e., most of 
the Mg lb lines), where the Gaussian profile was inappropri- 
ate. In a few cases the adopted EW was the average from Voigt 
and Gaussian fits. The Voigt fit is much more sensitive than 
the Gaussian fit to the continuum location. Any systematic er- 
ror in the continuum fitting would introduce a trend in the Mg 
lb abundances versus line-strength, and thus versus metallic- 
ity (the highest-metallicity stars showing the stronger lines). 
However, for these highest-metallicity stars we always have a 
weaker Mg line to check the derived abundances, as discussed 
below. 

To estimate the uncertainties on our measured EWs, we cal- 
culated the standard deviation for each of the 160 pairs of mul- 
tiple measurements, and obtained a mean standard deviation a 
-1.6 mA. This is in very good agreement with the 2 mA sensi- 
tivity suggested by RNB on the basis of the spectral S/N ratios. 
For the subset of stronger lines with EW > 70 mA , we found a 
mean standard deviation o^£vi'>70mA - 2.5 mA. The strong Mg 
lb lines have tf = 4 mA, corresponding to a 4% uncertainty in 
the derived EWs. 

3.3. Refining tlie line measurements 

We filter our line measurements by constraining the line central 
wavelength, measured equivalent width, and velocity width, 
FWHMv, corresponding to the Gaussian FWHM. We require 
the centre of the fitted profile to be within 0.05 A of the central 
wavelength observed in the solar spectrum (Moore et fl/.'1966't. 
This allows for an uncertainty in the adopted rest-frame correc- 
tion, and for a small deviation of the centre of the fit from the 
core of the line (e.g., as can occur when the fit is dominated 
by the wings of the profile). We then require a minimum EW 
of 7.0 mA, which is 4.5 times the mean standard deviation ob- 
tained for multiple observations, a". An analysis based on the 
more conservative criteria of requiring EW larger than 10.0 mA 
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shows an increased scatter in our final abundances; although 
the higher threshold leads to more precise measurements of in- 
dividual lines, the number of accepted lines decreases, in par- 
ticular for the most metal-poor stars. The lines were also fil- 
tered by requiring 7.0 < FWHMy < 14.0 km s"' in order to 
avoid blends, or lines so distorted by noise to appear either nar- 
rower than the instrumental profile or unphysically wide. This 
selection should lead to higher internal consistency of the mea- 
surements. The final EWs are reported in Tabled For each 
line we list the chemical species (el), central wavelength (/i), 
adopted ^/-value and coiTesponding reference (ref), and the 
equivalent width measured in each program star. 



3.4. Model Atmosphere 

We adopt the grid of ID, LTE model atmospheres of Kurucz 
(1989, private communication), interpolating within the orig- 
inal models to our required stellar parameters. As previously 
noted (Ryan et al. I1996> . abundances derived with Kurucz's 
1989 models are in good agreement with abundances derived 
with the model atmospheres of Bell et al. J1981> . but newer 
Kurucz (1993'i models have a different temperature structure 
and lead to abundances higher than the previous models. Ryan 
et al. (I1996t found that temperatures are up to 200 K higher 
in the line-forming region of the Kurucz ( 1993 1 models and 
result in -i-O.l dex in their Fe I abundances for turn-off stars. 
We find a +0.1 dex difference for lines of EW ~ 10 mA, ris- 
ing to +0.15 dex at ~ 40 mA and to +0.2 dex at ~ 150 mA, 
the wings of stronger lines being more sensitive to these re- 
gions. The two strong Mg lb lines are thus more sensitive to 
the change in the temperature structure of the model than the 
bulk of Fe lines. This would imply a typical -0.05 dex differ- 
ence in the [Mg/Fe] ratio if we derived the abundances with the 
Kurucz ( 1993 1 models and only the b lines of Mg I. However, 
the use of weak Mg lines allows us also to use estimates for 
the abundance from lines in the same line-forming region as 
the Fe lines, and thus obtain a [Mg/Fe] value that is less depen- 
dent on the temperature structure of the model. The average of 
Mg abundances from weak and strong lines reduces the depen- 
dence of our results on the adopted model. Moreover, for the 
lowest-metallicity stars (where weak Mg lines are not measur- 
able), the Mg lb lines are also weaker and thus less dependent 
on the temperature structure. We estimate that if we were to 
adopt Kurucz ( 1993) models, the standard deviation about the 
mean [Mg/Fe] would increase by less than 0.005 dex. We can 
thus safely neglect the model dependence in our discussion of 
the observed scatter. 



3.5. Stellar Parameters 

In principle, we can derive the effective temperature, Teff , in an 
iterative fashion by computing abundances for several models, 
and requiring no dependence of the Fe abundance on excita- 
tion potential (EP). The surface gravity log(g) can be set by 
requiring the Fe abundance derived from Fe II lines to match 
that from Fe I lines. Finally, the microturbulent velocity can 
be derived by requiring no trend of the abundances versus EW. 
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Fig. 1. The position in the diagram dereddened ci vs. b - y is 
shown as proof of the similar evolutionary stage of the stars. 
Filled squares represent stars with measurements of Fe II lines 
and thus spectroscopic estimate of log(g). Empty squares are 
used for stars with no Fe II measurements. 



However, these spectroscopic estimates for the stellar parame- 
ters rely on the measurement of at least some lines of both Fe I 
and Fe II. The higher the number of lines, the higher the preci- 
sion of the estimate. As discussed before, the very metal-poor 
stars show just four or five reliable lines in the stellar spectra; 
for about half of the sample we cannot measure any Fe II lines 
at all. Hence, we are not able to derive spectroscopic estimates 
of the parameters for all the stars and furthermore, when we do 
have estimates, not all have the same precision. 

Photometry can assist us in the task of determining the stel- 
lar parameters, having the excellent characteristic of not bias- 
ing our precision at lowest metallicity (as is the case for spec- 
troscopic estimates, which rely on a decreasing number of lines 
at lowest metallicity). As shown in Fig.[n the stars are in a re- 
gion of the C] vs. b - y diagram that led RNB to identify them 
as main-sequence turnoff stars. This determination is done by 
comparing the observations with evolutionary tracks, on the ba- 
sis of the b-y temperature sensitive index and the gravity sensi- 
tivity of the ci index. We go a step further, and obtain estimates 
for the surface gravity by calibrating the ci index with stars 
whose spectroscopic gravities we are able to derive. Because 
of the small number of Fe I lines and the absence of Fe II lines 
for the lowest-metallicity stars, we adopt photometric calibra- 
tions both for Teft and log(g), based on the homogeneity of the 
sample, rather than using solely spectroscopic values. We dis- 
cuss the details of our method below. 

Effective Temperature - 

The effective temperature is adopted from RNB. Their pho- 
tometric estimates were based on six indices (B-V, V-R, R-I, 
P, HP2 and b - y), calibrated with the theoretical calculations 



Table 1. Equivalent widths (mA) for program stars. 
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ref 


BD-13"3442 
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BBD+3"740 
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111 


Mg I 


5528.40 
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Cal 
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17 
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Cal 
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12 






















Cal 
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-0.31 


13 


16 
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25 






15 


12 
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13 


21 


22 


35 


35 


16 


42 


21 


15 




13 


Cal 
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12 






















Cal 
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12 
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10 
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5206.04 


0.02 


10 




17 


25 


21 




27 








14 


Crl 


5208.43 


0.16 


10 




21 


35 


34 




36 










CrI 


5409.80 


-0.72 


10 
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11 
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-2.87 


11,2,7 












10 










Fel 


5166.28 


^.20 


5 






















Fel 


5171.61 


-1.76 


11,6 


11 


14 


16 


28 


15 


26 


12 


16 




8 


Fel 


5191.47 


-0.55 


11 






12 


12 














Fe I 


5192.35 


-0.42 


11 




9 


15 






21 










Fe I 


5194.95 


-2.06 


11,6 
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10 




14 










Fe I 


5198.72 


-2.1 1 
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Fe I 


5269.55 


-1.33 


1 1 


57 


54 


66 


73 


61 




60 


54 


36 


42 


Fe I 


5281.80 


-0.83 


11 












1 1 










Fel 


5283.63 


-0.48 


11,2 


7 




11 


10 




14 




7 






Fel 


5302.31 


-0.72 


2 






















Fel 


5324.19 


-0.10 


2 




13 


20 


24 




27 


15 
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5328.05 


-1.47 


11, 7 


48 


44 


58 


68 
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48 


43 
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-1.85 


11 
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16 


15 




23 
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-1.95 
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22 
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34 


37 


24 
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23 


11 


14 


Fel 


5405.79 


-1.86 


11,2,7 


27 


22 


34 


40 


28 




25 


23 


9 


15 


Fel 


5410.92 


0.40 


11 
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Fel 


5415.21 


0.64 


11 


9 




10 


13 




16 










Fel 


5429.71 


-1.88 


11 


28 


20 


35 


44 


32 


46 


28 


22 


11 


18 


Fel 


5434.53 


-2.12 


11,2,7 


17 


13 


24 


25 


17 


30 


16 


13 




10 


Fel 


5446.92 


-1.88 


11 




18 


34 


38 


28 


41 


23 


24 


13 


11 


Fel 


5569.63 


-0.49 


1 












10 










Fel 


5572.85 


-0.28 


2 






9 






16 










Pel 


6136.62 


-1.41 


11 












11 










Fel 


6191.57 


-1.51 


11 












9 










Fel 


6678.00 


-1.42 


11 












7 










Fe II 


4923.93 


-1.32 


9 




38 


46 








36 


36 






Fen 


5197.58 


-2.10 


9 












15 










Fen 


5276.00 


-1.95 


9 








15 






10 


10 







1 Bard & Kock 1994, 2 Bard, Kock & Kock 1991, 3 Blackwell et al. 1982a, 4 Blackwell et al. 1982b, 5 Blackwell et al. 1982c, 6 Blackwell et al. 1980, 7 Blackwell et al. 1979, 8 Fuhrmann et 

al. 1995, 9 KroU & Kock 1987, 10 Martin et al. 1988, 11 O'Brian et al 1991, 12 Smith & Raggett 1981, 13 Wiese & Martin 1980. 
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of Bell & Oke (I1986> and the empirical relations of Magain 
J1987> . The photometric indices were de-reddened using the 
reddening maps of Lucke 1978 and Burstein & Heiles T5^), 
with Johnson photometric distances from Carney et al. 81994 ) 
and Ryan et al. (1989 1, and the comparison of b - y with 
the reddening-free index /3 for Stromgren indices. A correc- 
tion was applied to account for a 0.020 mag difference be- 
tween the colour excess estimate derived from Stromgren in- 
dices and that inferred from reddening maps, where we expect 
E{b - y)=0.7E(B- V). The uncertainties on the photometric in- 
dices were assumed to be on average less than +0.01 magni- 
tudes (see RNB and references therein for details). The final 
Teff estimates are shown in Table |2] They exhibit a very low 
relative uncertainty, on the of order ±40 K, although larger sys- 
tematic errors could affect the zero point by up to ~ 100 to 200 
K. 

Later spectroscopic reanalysis of some stars in our sam- 
ple have suggested higher temperatures (see e.g., Nissen et al. 
1200 It Ford et al. .2002) . We note the problem of systematic 
differences between photometric and spectroscopic estimates 
of Teff, the former being typically expected to be 100 K hot- 
ter according to Alonso et al. il999i) . Johnson (,2002) finds up 
to 150K difference between Tefj estimates derived with spec- 
troscopy or with photometry. We therefore determine a spec- 
troscopic estimate of the T^ff for the stars by requiring no de- 
pendence of the abundance derived from Fe I lines versus the 
excitation potential (BP) of the line. For most of the stars, there 
are too few Fe I lines to have a high confidence in the deter- 
mination of spectroscopic Tspec; we work on the subset of 15 
stars with more than 10 lines. Our estimates for the effective 
temperature Tspec are on average 80K lower than the tempera- 
tures derived by RNB. This is just at the limit of our average Icr 
statistical error of 80K in the derivation of Tspec- However, in 
most cases we have only two or three Fe I lines with high EP, 
hence the procedure is heavily dependent on these few lines. 
We thus rely instead on the photometric estimates from RNB, 
and assess the impact of a possible systematic error up to 150K, 
as suggested by the largest differences between RNB estimates 
and our spectroscopic determinations. We assume a 40K ran- 
dom error on the temperature from RNB. 

Surface gravity and microturbulent velocity - 

We adopt an iterative procedure to determine both the sur- 
face gravity and the microturbulent velocity of each star. We 
make calculations with log(g)=4.0 and ^=1.0-1.5 km s"' as 
our first stage. With these results we calculate the best estimate 
for ^ by requiring no trend of the abundances of Fe I with the 
EW of the line the abundance was derived from. Due to the lack 
of lines for the most metal-poor stars, we adopt a range for ^ of 
1.0-1.5 s ' to compensate for the poor spectroscopic determi- 
nation. We perform calculations with the assumed log(g)+0.2 
dex and the calculated ^, searching for the value of log(^) that 
would let the Fe abundances derived from Fe I match the ones 
from Fe II lines. We then iterate this procedure to obtain new 
estimates of ^ and log(g), finding a convergence of the results 
after a few iterations. 

Only one or two Fe II lines are available for each log(g) 
determination. The standard error in the difference [Fe/H]i- 
[Fe/H]u corresponds to an error criog(g)p^ ~ 0.12 dex. The in- 



Table 2. Model parameters 



star 


-^RNB 






l02f£)«npc 

^"to Vi^ /spec 


l0g(?)nhnf 








(±0.2) 




(±0.2) 


BD-13°3442 


6210 


30 


1.5 




3.6 


BD+r2341p 


6260 


40 


1.5 


3.8 


3.8 


BD+20°2030 


6200 


40 


1.6 


3.9 


3.9 


BD-l-24°1676 


6170 


30 


1.8 


3.7 


3.7 


BD-l-26°2621 


6150 


40 


1.6 




4 


BD-i-26°3578 


6150 


40 


1.5 


3.4 


3.7 


BD-l-3°740 


6240 


40 


1.8 


3.7 


3.7 


BD-l-9°2190 


6250 


30 


1.3 


3.8 


3.6 


CD-24° 17504 


6070 


30 


1.8 




4.2 


CD-33°1173 


6250 


20 


1.8 




3.7 


CD-35° 14894 


6060 


20 


1.3 


4.0 


4.1 


CD~71°1234 


6190 


30 


1.5 


4.1 


4.1 


CS22943-045 


6140 


40 


1.5 


3.7 


3.9 


G 126-52 


6210 


40 


1.8 


3.9 


3.8 


G4-37 


6050 


40 


1.8 




4.1 


G64-12 


6220 


30 


1.7 




3.9 


G64-37 


6240 


30 




3.9 


3.9 


HD 74000 


6040 


30 


1.3 


4.2 


4.1 


HD84937 


6160 


30 


1.8 


3.9 


3.7 


LP635-14 


6270 


30 


1.7 


4.4 


3.7 


LP651^ 


6240 


30 


1 




3.9 


LPS 15-43 


6340 


30 


1 




3.6 


LP831-70 


6050 


20 


1 




4.3 



ferred log(g) value is also affected by the adopted effective 
temperature. Changes in T^ff of 40 K cause a variation in the 
inferred \og{g) of cTia^fgyj. ~ 0.02 dex. Values generally adopted 
for ^ in the case of turn-off stars are in the range 1.0 - 1.5 
km s"'. A conservative assumption of a variation of ^, over 
the range A^=± 0.2 km s produces a variation in the de- 
rived \og{g) of ± 0.06 dex on average. On the other hand, a 
variation of log(^) of ± 0.06 produces a variation in ^ of less 
than 0.03 km s Thus the determination of ^ is quite insen- 
sitive to the adopted \og{g). We expect this 0.06 dex uncer- 
tainty on \og{g) to arise only in the worse cases when we can- 
not constrain ^; in most cases this iterative procedure produces 
an uncertainty in the derived \og{g) of criog((,)p^^^^j^^^ ~ 0.04 dex. 
We thus assume a random error on the spectroscopic \og(g) of 

o"ioe(e) = -./o"? ,^ + o"? , , + ai , . =0.13. For 4^ we 
adopt the more conservative uncertainty of 0.2 km s ' . 

These spectroscopic estimates of \og{g) are used to cali- 
brate logig) versus the de-reddened ci index using a linear 
fit (see Fig. |2}. The ci index was de-reddened adopting the 
colour excess law E{c\) = 0.20E{B - V) and ci,o - ci - E(ci) 
(Golav.1974). The calibration is then used to calculate 'pho- 
tometric' estimates of the logig) for each star (Table |2j). The 
standard deviation about the fit for the calibration of log(§) ver- 
sus ci is 0.15 dex, excluding (by a 2cr clipping procedure) the 
log(g)=4.4 derived spectroscopically from just one Fe II line 
for LP635-14. The error on the photometric estimates of log(§) 
is adopted to be 0.2 dex, as the quadratic sum of the uncertainty 
on the spectroscopic log(g) (0.13) and of the calibration versus 
ci (0.15), the influence of the uncertainty on the ci value be- 
ing negligible. We then derive the abundances with the adopted 
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Fig. 2. Calibration of the spectroscopic log(^) vs the gravity 
sensitive dereddened ci photometric index. The empty square 
represents the star LP635-14 that was rejected from the fit by a 
2cr clipping procedure. 

photometric \og{g) for the entire sample of stars. For the subset 
of stars for which we could calculate the spectroscopic \og{g), 
we perform a further analysis with that value of \og{g), but find 
no significant difference in the overall trend or scatter of the 
results. 

We investigate a possible dependence of the ci index on 
temperature as well as gravity, but find no evidence of a corre- 
lation between the residuals of the gravity calibration and the 
temperature scale. This indicates that the gravity sensitivity of 
the Ci index is far more important than its temperature sensi- 
tivity over the narrow temperature range of our stars, and the 
latter does not aff'ect our calibration. 

Initial chemical composition - 

We also explored the possible sensitivity of the derived 
abundances on the initial chemical composition adopted for 
the calculations. The chemical composition is introduced in the 
model atmosphere by scaling the solar composition with the 
metallicity of the star However, because of the very low metal- 
licity of our stars, no significant change in the derived abun- 
dances is produced by adopting models with slightly different 
chemical composition. 

4. Chemical abundances 

We now derive surface abundances from lines of Mg I and Fe 
I for our entire sample of stars. For most of the stars we are 
also able to derive surface abundances for the elements Ca, Ti, 
and Cr, depending on the presence of reliable spectral lines in 
our wavelength coverage. The results are listed in Tables|3]and 
13 along with the number of lines used, as well as the standard 



deviation, cr, of the abundances derived from different lines. 
The total uncertainty of each element is described below. 

4.1. Dependence of abundances on model 
atmosphere parameters 

The dependence of our derived abundances on stellar parame- 
ters varies with the strength of the lines used. Because of the 
homogeneity of our sample, we study this dependence for the 
whole sample rather than for single stars. 

For Fe I, the number of weak lines is far higher than the 
number of stronger ones. This reduces the sensitivity to strong 
lines, hence there will be no large differences of sensitivity be- 
tween higher and lower metallicity stars. Fe II lines are never 
stronger than 60 mA , so again the averaged Fe II sensitivity 
is a good representative for all the stars. We derive abundances 
for a grid of models in which we vary independently one stellar 
parameter at a time by its standard error, as discussed in Sec. 
13.51 and report the average results for Fe lines in Table|5] 

In the case of Mg, the model parameter sensitivity is much 
more significant in the Mg lb lines than in the weak Mg line. 
Besides a dependence on temperature similar to that of Fe I 
lines, Mg lb lines show an extremely high sensitivity to surface 
gravity because of the stronger wings. We show in Table|5]the 
average results for the weak Mg I line at 5528 A and for the 
average Mg lb lines. The dependence of the [Mg/Fe] ratio is 
averaged on the whole sample of stars, and thus weighted on 
both the weak line and the two strong lines. 

Estimating uncertainties by varying just one parameter at a 
time is a useful exercise, but it does not reflect the real anal- 
ysis procedure. In order to estimate more realistic uncertain- 
ties, we performed a Monte Carlo simulation reproducing the 
procedure adopted for deriving \og{g) and ^. We consider two 
candidate stars, BD-i-3°740 (10 Fe I fines, [Fe/H] = -2.84) and 
BD-i-26°3578 (21 Fe I fines, [Fe/H] = -2.49) as representative 
of our sample. For each star we create a set of 25 virtual "ob- 
servations" by adding numerical noise to the measured EWs. 
For each line, the distribution of the 25 randomly generated 
eiTors added to the EW is a Gaussian with cr = 1.5 mA (i.e., 
our average inferred Icr uncertainty on the EW). By deriving 
the surface gravity and microturbulent velocity for the virtual 
observations with the iterative method used for the real data, 
we were able to estimate the dependence of our abundances 
and derived parameters on observational errors, by means of 
the standard deviation of the results for the 25 simulations. In 
Table |6l we report the results of this "end-to-end" test of our 
procedures. 

From this test we find that the sensitivity of our results 
is lower than a quadratic sum of the single parameter depen- 
dences estimated one at a time. This is due to the fact that 
the errors in the parameters derived with the iterative proce- 
dure are not independent, so the parameters partially compen- 
sate for one another Moreover, the abundances are averaged 
on both weak and stronger lines, decreasing the sensitivity that 
strong lines would show alone. It is comforting to notice that 
the standard deviation of the ^ and log(g) values are in the or- 
der of 0.2 km s ' and 0.2 dex respectively, as expected from 
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Table 3. Fe and Mg chemical abundances. 
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star 


[Fe/H]i 


npel 




[Fe/H]„ 


npeii 


Fell 


[Mg/H] 






[Mg/Fe] 


BD-13°3442 


-2.83 


9 


0.02 


— 








-2.47 


3 


0.03 


0.36 


BD+r2341p 


-2.88 


11 


0.02 


-2.85 


1 





-2.51 


3 


0.01 


0.37 


BD+20°2030 


-2.68 


16 


0.01 


-2.68 


1 





-2.39 


3 


0.06 


0.29 


BD+24°1676 


-2.60 


15 


0.02 


-2.60 


1 





-2.28 


3 


0.01 


0.32 


BD+26°2621 


-2.86 


7 


0.01 


— 








-2.59 


3 


0.01 


0.27 


BD+26°3578 


-2.49 


21 


0.01 


-2.40 


1 





-2.16 


2 


0.01 


0.33 


BD+3°740 


-2.84 


10 


0.01 


-2.87 


2 


0.09 


-2.52 


3 


0.03 


0.32 


BD+9°2190 


-2.83 


10 


0.02 


-2.87 


2 


0.07 


-2.61 


3 


0.07 


0.22 


CD-24° 17504 


-3.45 


5 


0.03 


— 





— 


-3.17 


3 


0.10 


0.28 


CD-33°1173 


-3.10 


8 


0.03 


— 





— 


-2.81 


3 


0.04 


0.29 


CD-35° 14894 


-2.55 


20 


0.01 


-2.53 


2 


— 


-2.31 


2 


0.01 


0.24 


CD-71°1234 


-2.55 


22 


0.01 


-2.56 


3 


0.01 


-2.35 


3 


0.03 


0.20 


CS22943-045 


-2.50 


22 


0.02 


-2.43 


2 


0.08 


-2.13 


3 


0.02 


0.37 


G126-52 


-2.44 


20 


0.01 


-2.49 


1 


— 


-2.17 


2 


0.05 


0.27 


G4-37 


-2.75 


12 


0.02 


— 





— 


-2.48 


2 


0.00 


0.27 


G64-12 


-3.35 


4 


0.04 


— 





— 


-3.08 


3 


0.05 


0.27 


G64-37 


-3.15 


8 


0.02 


-3.13 


1 




-2.92 


3 


0.04 


0.23 


HD 74000 


-2.17 


28 


0.01 


-2.24 


1 




-1.95 


3' 


0.07 


0.22 


HD84937 


-2.36 


23 


0.01 


-2.44 


3 


0.01 


-2.01 


3 


0.02 


0.35 


LP635-14 


-2.53 


23 


0.02 


-2.80 


1 




-2.19 


3 


0.02 


0.34 


LP651-4 


-2.69 


14 


0.03 







0.04 


-2.27 


3 


0.04 


0.42 


LP8 15-43 


-2.81 


7 


0.01 









-2.42 


2 


0.09 


0.39 


LP831-70 


-3.06 


9 


0.03 









-2.69 


3 


0.06 


0.37 



' : In HD74000 we measured the Mg weak line at 571 1 mA instead of the line at 5528 mA. 



the procedure used in Sec. 13. 51 More interestingly, the resulting 
dependence of the abundances on the EW uncertainties and, in 
turn, on these two imperfectly determined parameters, is very 
small, of the order 0.02 dex for [Fe/H] and 0.04-0.06 dex for 
[Mg/H] . The real impact of the parameter determination is thus 
much smaller than Table |5] suggested. This procedure is also 
useful for evaluating how these errors cancel out in the ratio 
of the two elements, since we are interested in determining the 
[Mg/Fe] ratio; the uncertainty varies between 0.03 dex and 0.06 
dex for stars of higher and lower metallicity. We re-emphasize 
that we did not derive T^s spectroscopically in this process, 
because of the lack of lines of high BP. Consequently, an un- 
certainty from the adopted Teff must still be included, which we 
take from Table|5] However, this uncertainty almost cancels out 
for the [Mg/Fe] ratio. 

4.2. Non-LTE effects on the abundances 

We derive chemical abundances assuming local thermodynam- 
ical equilibrium (LTE). This assumption may be unrealistic, 
and one should be careful that lines may form in conditions 
that depart from LTE. However, the current efforts to remove 
the LTE assumption yield conflicting results, as shown by the 
results of Thevenin & Idiart ( 1999 1 and Gratton et al. (1999 1 
for Fe, and Zhao et al. (IT998i for Mg I. 

Non-LTE (NLTE) calculations from Thevenin & Idiart 
J199 9l show that iron abundances derived from Fe I lines 
should be increased by 0.3 dex, compared to the LTE esti- 
mates, for [Fe/H] in our metallicity range. On the other hand, 
Mg abundances should be reduced by 0.04 dex when derived 



from the Mg I line at 5528 A and by 0.03 dex for the two Mg 
lb lines (Zhao et al. ,1998j. This would lead to a systematic 
overestimate of [Mg/Fe] of 0.33 dex in LTE (with a 0.04 dex 
diff'erence in the NLTE corrections between the lowest and the 
highest metallicity of our stellar sample). However, since Fe is 
mostly present in its ionized state, Fe abundances derived from 
Fe II lines are insensitive to overionization, and a higher value 
of \og{g)NLTE (by 0.3-0.4 dex following Thevenin & Idiart 
(1999)) would be required to obtain ionization balance. This 
would further increase the overestimate of [Mg/Fe] to 0.45 dex 
in LTE. 



It is important to note that Gratton et al. (•1999^ found dif- 
ferent results when calculating NLTE corrections to Fe and Mg 
abundances. They suggest corrections that are very small, if not 
negligible, for main-sequence and RGB stars. This is in con- 
flict with the calculations of Thevenin & Idiart U999J . leaving 
NLTE corrections still very uncertain. 



The assumption of LTE is thus presumably introducing er- 
rors in our abundances but, because of the homogeneity of our 
sample and the small variation of NLTE corrections on the 
metallicity range considered, we assume any effect on [Mg/Fe] 
is buried in our random uncertainties. Furthermore, because of 
the present uncertainties on the effects of NLTE, and in order 
to have results comparable with most observational studies, we 
adopt the LTE estimates in our discussion, warning the reader 
that our [Mg/Fe] ratio may be overestimated. 
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Table 4. Ca, Ti and Cr chemical abundances. Ti II abundances are derived only in 3 stars. 
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star 


[Ca/H] 




0"Ca 


[Ti/H] 


nxi 




[Cr/H] 


nc, 


o"c, 


BD-13°3442 


-2.4 


3 


0.03 


— 





— 


— 





— 


BD+r2341p 


-2.37 


3 


0.04 


-2.25 


1 


— 


-2.84 


2 


0.01 


BD+20°2030 


-2.23 


3 


0.04 


-2.16 


1 


— 


-2.62 


2 


0.04 


BD+24°1676 


-2.21 


3 


0.03 


-2.07' 


1 


— 


-2.72 


2 


0.08 


BD+26°2621 


-2.67 


1 


— 


— 





— 


— 





— 


BD+26°3578 


-2.05 


2 


0.02 


-2.10 


1 


— 


-2.62 


2 


0.04 


BD+3°740 


-2.43 


2 


0.01 


— 





— 


— 





— 


BD+9°2190 


-2.51 


3 


0.06 


— 





— 


— 





— 


CD-24° 17504 


— 





— 


— 





— 


— 





— 


CD-33°1173 


-2.69 


1 


— 


— 





— 


-2.96 


1 


0.0 


CD-35° 14894 


-2.22 


1 


— 


— 





— 


-2.69 


2 


0.07 


CD-7 1-1234 


-2.23 


3 


0.02 


— 





— 


-2.51 


3 


0.18 


CS22943-045 


-1.99 


5 


0.05 


-2.08^ 


1 


— 


-2.57 


2 


0.04 


G 126-52 


-2.09 


2 


0.00 


-2.07 


1 


— 


-2.54 


1 


— 


G4-37 


-2.39 


1 


— 


-2.36 


1 


— 


-2.81 


2 


0.04 


G64-12 


-2.87 


1 


— 


— 





— 


— 





— 


G64-37 


-2.86 


1 


















HD 74000 


-1.93 


4 


0.05 


-1.91 


2 


0.03 


-2.24 


3 


0.08 


HD84937 


-1.94 


4 


0.04 


-1.97' 


1 




-2.44 


2 


0.04 


LP635-14 


-2.22 


3 


0.02 


-1.78 


1 




-2.63 


2 


0.04 


LP651-4 


-2.43 


2 


0.07 
















LP8 15-43 


-2.47 


1 


















LP831-70 









-2.46 


1 




-3.20 


1 




[Ti/H]n = 


-2.19,2 


[Ti/H]„ 


= -2.09 and ' 


[Ti/H]„ 


= -2.03. 







Table 5. Dependence of abundances on model parameters 



parameter 


variation 


A[Fe/H], 


A[Fe/H]„ 


A[Mg/H]j528 


A[Mg/H]ib 


A[Mg/Fe] 


Teff {random) 


±40K 


±0.03 


±0.00 


±0.02 


±0.04 


±0.00 


T^ff (systematic) 


±150K 


±0.13 


±0.02 


±0.07 


±0.14 


±0.01 




±0.2 


+0.01 


±0.00 


+0.01 


+0.03 


+0.02 


log(g) 


±0.2 


+0.01 


±0.03 


+0.02 


+0.10 


+0.06 



Table 6. Standard deviations for the results of the Monte Carlo test 



star 


[Fe/H] 


cTEw/mA 




criog(g) 


O" [Fc/HJ, 


cr [Fc/H]„ 


cr [Mg/H] 


cr [Mg/FcJ 


cr [Ca/Fe] 


BD+3°740 


-2.84 


1.5 


0.20 


0.19 


0.02 


0.03 


0.06 


0.06 


0.02 


BD+26°3578 


-2.49 


1.5 


0.17 


0.19 


0.02 


0.03 


0.04 


0.03 


0.02 



4.3. Abundance error estimates 



4.3.1 . Random errors 



In Sec. l3.5l we estimated the uncertainties on the adopted stel- 
lar parameters. In Table |5] we showed the dependence of the 
abundances on these parameters taken one at time; in Table |6l 
the results of the Monte Carlo test showed the combined de- 
pendence on iteratively determined log(g) and f values to have 
a smaller impact on the abundances than the quadratic sum of 
the two errors taken separately. Moreover, the Monte Carlo test 
was performed propagating numerically the uncertainty on the 
EW, thus including the contribution of the latter in the result. 
We can now estimate the errors on our final abundances. We 
treat separately the random errors that affect the scatter of our 
results, and the systematic errors that affect the absolute values 
(and thus the average [Mg/Fe]). 



As shown in Table |6l the uncertainty on the EW propagates 
through the determination of the parameters log(g) and ^, and 
affects the resulting [Fe/H] by 0.02 dex, [Mg/H] by up to 0.06 
dex, and [Mg/Fe] by up to 0.06 dex. We could adopt these val- 
ues as the uncertainties on the abundances due to the combina- 
tion of the three factors, lcrj^{crEw,criog(g),cr^). However, we 
make use of log(g) values from the photometric calibration, 
hence the Monte Carlo test does not reflect the entire proce- 
dure. We are obliged to estimate our errors as the quadratic 
sum of the uncertainty due to single parameter variations, and 
consider the results of the Monte Carlo test as a conservative 
estimate of the error due to uncertainty on line measurement, 

O-AicTEw)- 
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Adding the random uncertainty on temperature, ^ and 
log(g) (Table |5} increases the uncertainty on the abundances 
to 0.04 dex for [Fe/H], 0.10 dex for [Mg/H]ib and 0.08 dex for 
[Mg/Fe]n,. Stars where the weak Mg 5528 A line was not mea- 
sured have a slightly higher uncertainty on Mg, by about 0.01 
dex. 

Log (gf) values have been chosen on the basis of an exper- 
imental determination with uncertainties <25%. This implies 
that we have to account for a coiTesponding < 0. 1 dex uncer- 
tainty in the abundances derived from single lines. Having a 
number of lines will reduce the uncertainty in the average by 
a factor 1/ V^V. For Fe abundances, where the number of lines 
is statistically significant, we assume the standard eiTor of the 
Fe abundance averaged on all the lines to be a good estimate 
of the total random error. This is up to 0.03 dex. For Mg, the 
two or three lines give a random error due to ^/-values of ~ 
0.06 dex. However, in the case of Mg, we treat these errors as 
systematic, as described in the next section. Adding these un- 
certainties quadratically to the previous estimates produces a 
total random error of 0.05 dex for [Fe/H], 0.10 dex for [Mg/H] 
and 0.08 dex for [Mg/Fe] . 

A final remark is relevant on the uncertainties derived by 
RNB on their temperature estimates. Since their 40K random 
uncertainty is very low, and our [Fe/H] abundances are very 
sensitive to temperature, we explored the effect of adopting an 
artificially more conservative cr = 80K error for this study. This 
would increase the random uncertainty on [Fe/H] to 0.07 dex 
and that on [Mg/H] to 0.11 dex. The [Mg/Fe] error remains 
unchanged at 0.08 dex. 

4.3.2. Systematic errors 

We estimate the effect of systematic errors in our model param- 
eters as follows. We consider the effect of an overestimate of 
150 K as the systematic eiTor on the temperature, based on the 
previous discussion of the photometric estimate adopted in this 
study — see Sec. 13. 51 The variation on [Mg/Fe] for a systematic 
change of -150 K alone would be less than -0.01 dex, because 
of an equivalent decrease of ~ 0.14 dex in both [Fe/H] and 
[Mg/H] . However, the procedure used to determine the spectro- 
scopic log(g) would partially compensate for the Tefj change, 
resulting in a decrease in log(g) of about 0.3 dex, leading to an 
overall increase of 0.04 dex for [Mg/H]. The overall variation 
on [Fe/H] is 0.12 dex, because of the much weaker sensitiv- 
ity to changes in log(g). This produces a 0.08 dex change in 
[Mg/Fe], because of the greater sensitivity of Mg lb lines to 
log(§). 

The uncertainty on ^/-values can be a large contribution to 
our uncertainty on the abundances, and requires further com- 
ment. For Fe, we introduced the uncertainties on the ^/-values 
in the random error. However, in the case of the Mg lb lines, 
the uncertainty on the ^/-values should be treated like a sys- 
tematic effect, since the two Mg lb lines are averaged only with 
the weak Mg line at 5528 mA. Because of the similar EPs and 
strength of the two Mg lb lines, eiTors in the Mg ^/-values 
could cause the same variation in the line abundances of both, 
thus shifting our Mg abundances by the same amount. This will 



RNB 1 A estimate 
literature 




-3.5 - 



-3.5 -3 -2.5 -2 

[Fe/H].„.„,, 

Fig. 3. A comparison between the Fe estimates of this study 
and the 1 A (filled squares) and literature (empty squares) esti- 
mates respectively. Overplotted are the two linear fits. 



affect the mean value of the Mg abundance of a star but not the 
line-to-line scatter about the mean. We prefer to include the un- 
certainty in the Mg ^/-values separately from the random error 
derived in the previous subsection. If we do so, the random er- 
ror on [Mg/Fe] is 0.08 dex, while the systematic eiTor becomes 
the sum of the systematic error due to temperature (0.08 dex) 
and the one due to the ^/-values (0.06 dex), i.e. the total sys- 
tematic uncertainty on [Mg/Fe] is then 0. 10 dex. The total sys- 
tematic error for [Mg/H] is 0.07 dex, while [Fe/H] is affected 
systematically only by the temperature, i.e. by 0. 12 dex. 

4.4. Metallicity comparison 

Figure|3shows a comparison of the metallicity estimates of this 
study with the 1 A estimates of RNB (filled dots). The lA esti- 
mates were obtained by a calibration of the pseudo-equivalent 
width of the Ca II K line as a function of B-V colour to indi- 
cate the stellar metallicity [Fe/H] (Beers et al.'1999't. A linear 
fit ax + b gives a = 1.05 and b - 0.23, with an rms of 0.14 
dex, larger than our uncertainties on Fe (0.05), but consistent 
with the value of 0.15 dex claimed by RNB for the 1 A esti- 
mates. In addition, we show the trend of Fe literature values as 
discussed by RNB (empty squares). A linear fit gives a = 0.96 
and b = -0.16, with a smaller scatter about the fit, rms=0.11 
dex. This implies a slightly better agreement of our metallicity 
determination with previous Uterature determinations than with 
the 1 A estimates. Both, however, can be considered in good 
agreement with our estimates within the stated uncertainties. 
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5. Discussion 

5.1. Fe and Mg abundances 

The abundances of Fe and Mg are shown in Table [Mg/Fe] 
is within a range of 0.20 dex, between 0.22 and 0.42. Figure 
13 displays the run of [Mg/Fe] vs. [Fe/H]. The thick line is the 
result of a linear least squares fit and the two thin lines delimit 
a +lcr region, cr = 0.06 dex being the standard deviation of 
points about the fit. The equation of the linear fit is: 

[Mg/Fe] = (0.00 + 0.04) ■ [Fe/H] + (0.31 + 0.12) (1) 

The result of this linear fit suggests that the trend of 
[Mg/Fe] as a function of [Fe/H] is flat. We calculate a weighted 
mean for [Mg/Fe], the weights being the inverse of the square 
of the adopted uncertainties. The weighted average, standard 
eiTor and standard deviation are 

[Mg/Fe] = 0.30 + 0.01 cr = 0.06 (2) 

Systematic eiTors are not accounted for in this statement; 
they may affect the mean value by up to 0.10 dex, which we 
take to be the limit on the precision with which the average 
value of [Mg/Fe] is determined. 

The scatter is not affected by systematic errors, hence, if 
no real intrinsic scatter was present, we would expect the rms 
scatter (0.06 dex) to be no larger than the random uncertainty 
(0.08 dex), as is the case. Moreover, the standard deviation of 
our results is similar to that of the Monte Carlo simulation that 
was performed starting from the data of a single star (but which 
excluded random eiTors in Tetj). We thus infer that the disper- 
sion in the results is due to our analysis; any intrinsic scatter in 
[Mg/Fe] of the sample must be much smaller, <& 0.06 dex. This 
is our most important result, and is a considerable reduction on 
the Cayrel et al. standard deviation on [Mg/Fe], 0.13 dex. We 
consider the implications in Sec. 15.31 

The very low scatter confirms that our selection criteria and 
analysis produced a very homogeneous sample. A higher scat- 
ter could have been the signature of an intrinsic characteristic of 
the stellar sample or star-to-star differences introduced by the 
analysis. The latter can be difficult to estimate because of its de- 
pendence on the adopted model atmospheres. As an example of 
the difficulties in quantifying the real scatter, Johnson's (2002i 
study of giants has a dependence of [Mg/Fe] on both log(g) and 
T, [Mg/Fe] showing a very low scatter in the run versus each 
parameter but a high scatter versus [Fe/H]. A similar problem 
is buried in the results of Cayrel et al. (2003,) whose mean, 
[Mg/Fe] = +0.27 dex, is in reasonable agreement with our find- 
ings, but their standard deviation, 0.13 dex, has a component 
due to a slight dependence of their abundances on temperature 
and gravity. Our utilisation of turnoff stars covering only a nar- 
row range of temperatures and surface gravity has avoided this 
difficulty. In less homogeneous studies, different stellar evolu- 
tionary stages are considered in the same sample, thus includ- 
ing a star-to-star difference in the dependence on the adopted 
stellar atmospheres. Finally, the use of similar absorption lines 
for all the stars has limited the influences of uncertainties in gf- 
values to a systematic uncertainty on the average [Mg/Fe], not 
affecting the scatter. Further studies of larger, critically-selected 



samples will aid in determining the intrinsic properties of the 
Galaxy. 

5.2. Other elements 

The abundances derived for the other elements are shown in 
Tablel^based on up to five lines of Cal, three of Ti 1 and three 
of Cr 1. The trends of these abundances are shown in Fig|3 

For [Ca/Fe], we show a linear least square fits as a func- 
tion of [Fe/H]. A statistical test indicates that the slope is not 
significant. Since we have Ca estimates for all but two stars of 
our sample (CD-24° 17504 and LP831-70), we use [Ca/Fe] to 
verify the consistency of our ff-element abundances. The mean 
[Ca/Fe] is +0.37, with a slightly higher standard deviation than 
for [Mg/Fe], 0.09 dex. The estimated random error on [Ca/Fe] 
is slightly lower than the one on [Mg/Fe] - on average 0.07 
dex. (The Monte Carlo simulation, which excludes the temper- 
ature and the photometric calibration for the gravity, showed a 
much smaller error for [Ca/Fe], of order 0.02 dex, because of 
the smaller sensitivity of the weak lines used for Ca to changes 
in stellar parameters.) Consequently, we cannot rule out the ex- 
istence of an intrinsic scatter in the [Ca/Fe] values up to o" ~ 
0.06 dex. 

Previous authors have used an a-index averaging Mg 
and Ca chemical abundances (see e.g., Gratton et al. 2000 1. 
However, the Thielemann et al. (1996 1 SN 11 models show 
that Ca yields are only slightly dependent on progenitor mass, 
whereas Mg yields depend quite strognly on progenitor mass. 
This leads iGCE models to predict a lower scatter for Ca than 
for Mg (Argast et al. 2000 ). Moreover, Ca is synthesised deeper 
in the progenitor than Mg, so Ca yields are dependent also on 
the energy of the explosion, and hence not as independent as 
possible of SN 11 model uncertainties. We therefore concen- 
trate on [Mg/Fe]. 

We also tabulate abundances for Cr and Ti. Cr absolute 
abundances are, on average, similar to those of Fe. We con- 
firm the well-known overabundance of Ti relative to Fe, though 
with a suspicious dependence on [Fe/H] . However, the Ti abun- 
dances have been poorly determined due to the use of only one 
line for most of the stars, and hence we place little weight on 
this result. We do not investigate this element further 

5.3. Implications for Galactic chemical evolution 
5.3.1 . Restatement of the problem 

The absence of intrinsic scatter in the Mg abundances presents 
a real challenge for inhomogeneous Galactic chemical evo- 
lution models. Argast et al. C2QQ0 2002) and Tsujimoto & 
Shigeyama J1998> have shown that iGCE models predict a 
large scatter in [Mg/Fe] at [Fe/H]< -2.5 that reflects the de- 
pendence of nucleosynthesis on the progenitor mass of individ- 
ual SN II, at a time in the history of GCE when the ISM was 
poorly mixed. The Fran9ois et al. (2004") homogeneous GCE 
model can reproduce the abundance patterns observed in the 
most recent observations by requiring a large correction to SN 
II yields. One has to investigate why homogeneity occurs in a 
metallicity region expected to be inhomogeneous. 
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Fig. 4. Abundances are shown as [el/Fe] vs [Fe/H]. Empty squares are stars with no Fe II measurements. In the case of Mg and 
Ca, we show the results of a linear fit (thick line) and the +lcr region. This is not shown for Ti and Cr because of the small 
number of stars with Ti and Cr measurements. 



The Argast et al. model predicts that for [Fe/H]< -3.0, 
SN II pollute only locally, i.e., the ISM is mixed at a rate 
slower than successive SN II events, and individual SN yields 
are reflected by the local ISM. In the range -3.0 < [Fe/H] 
< -2.0, mixing of the ISM causes the first overlaps of pol- 
luted regions, though some regions that are unenriched remain. 
At [Fe/H]= -2.5, some regions of the ISM are enriched several 
times by supernova progenitors of different masses. One should 
thus observe stars showing abundance patterns characteristic 
of SN progenitors of different masses. By [Fe/H] > -2.0, the 
yields from successive polluting events are averaged over the 
initial mass function (IMF) and well mixed, bringing iGCE 
and 1-zone GCE models into agreement. Our metallicity range 
(-3.4 <[Fe/H]-2.2) covers the inhomogeneous part of the ISM 
history, where a range up to 1 dex is predicted for Mg abun- 
dances, compared to a range of 0.2 dex for our observations. 

In Figs.|5]and|6l we show the standard deviation of Argast's 
model calculated for 0.1 dex metallicity bins. The two fig- 



ures correspond to the models shown in Figs. 1 and 14 of 
Argast et al. (2002i, calculated adopting SN II yields from 
Thielemann et al. £l996l, and from the Argast et al. 'HI' em- 
pirical model, respectively. The yields from Thielemann et al. 
(19961 and Nomoto et al. (I1997> . adopted in the Argast et al. 
model, range from [Mg/Fe] =-0.83 for a progenitor of 1 3 M© to 
[Mg/Fe] = 1 .44 for a progenitor of 70 M©, with a monotonic re- 
lationship between the Mg mass yield and the progenitor mass. 
The HI model was obtained by modifying the SN II yields to 
match the Fe yields inferred from SN II observations, rather 
than using yields from a SN II model, and was further con- 
strained by the range of abundances previously observed in 
low-metallicity stars. This model shows a much lower scatter 
due to a smaller adopted range in SN II Mg/Fe yield ratios. 
Superimposed on these figures are shown our standard devia- 
tion for [Mg/Fe], and the predicted mean standard deviations 
of the models calculated over the same metallicity range. We 
also show the standard deviation of [Mg/Fe] from Cayrel et al. 
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-3.5 -3 -2.5 -2 -1.5 -1 

[Vc/H] 

Fig. 5. The standard deviation of the dispersion in [Mg/Fe] pre- 
dicted by the model of Argast et al. (.2002. their Fig. 1), in the 
case of SN II yields adopted from Thielemann et al. J1996t . 
The dotted, small dashed and the dashed lines show the stan- 
dard deviation over the metallicity range of our observations, 
-3.4 <[Fe/H]< -2.2, respectively for Argast et al. ( EUU^ . 
Cayrel et al. J2003I I and this work. 



J2003> . For the HI model, the predicted standard deviation in 
[Mg/Fe] is ~0.17 dex, rising towards lower metallicities, and 
well in excess of our limit of 0.06 dex. 

Theoretically (Arnett 1996), the assumption of instanta- 
neous mixing adopted by 1-zone GCE models is not correct 
at these low metallicities, hence inhomogeneity must be con- 
sidered in the models. If iGCE models are correct in relaxing 
the instantaneous mixing assumption, then the lack of scatter in 
the [Mg/Fe] values we have measured is all the more remark- 
able. To explain this result, one has to focus on SN II yields and 
other uncertainties in the models. 

5.3.2. A restricted mass range for SN II progenitors 

Chiappini et al. J 19991 1 found that, for metallicities 
-3.0 <[Fe/H]< -1.0, the trends in a-elements are ex- 
pected to show slopes. This is not seen in our flat trend 
at [Fe/H] < -2.2. However, this slope is not predicted in 
other models (e.g. Argast I2002t . The uncertainty on our 
measurement of the slope restricts the variation of the [Mg/Fe] 
slope to within 0.04 dex from flat over the metallicity range 
-3.0 <[Fe/H]< -2.0. 

The range of progenitor masses that give rise to 
[Mg/Fe] =0.3+0.1 (our observational value) is 18-20 Mq ac- 
cording to the yields of Thielemann et al. J1996> . However, 
even if we knew the exact yields of this progenitor mass range, 
the flat trend and the lack of intrinsic scatter raises the question: 
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Fig. 6. The standard deviation of the dispersion in [Mg/Fe] 
predicted by the model of Argast et al. (2002 their Fig. 14), in 
the case of their HI model. The dotted, small dashed and the 
dashed lines show the standard deviation over the metallicity 
range of our observations, -3.4 <[Fe/H]< -2.2, respectively 
for Argast et al. J2002> . Cayrel et al. ( I2003> and this work. 



Why, at any time in the early Galaxy, would the only contribu- 
tion to the ISM chemical composition of Mg come from this 
narrow range of progenitors? 

The upper mass limit that avoids complete collapse to 
a black hole is unknown, and while it is believed to be in 
the range of 30-50 Mq, there is some possibility it could be 
lower in low-metallicity stars. Maeder (1992 1 found that a limit 
around 20-25Mo was consistent with A7/AZ (the enrichment 
of helium relative to metals). More recently. Fryer J1999> con- 
firmed the 20 Mq limit for the initial mass for black hole for- 
mation. Fryer and Heger (2000 1 suggested that stellar wind and 
rotation could decrease the energy of the explosion and thus the 
lower mass limit for black hole formation. 

However, even if the upper mass limit for Mg enrichment 
of the ISM could be 20 Mq, there is no reason to think that pro- 
genitors as low as 10-12 Mq would be disabled. One may ques- 
tion our theoretical understanding of SN II models, or look for 
criteria that would favour SN II with this progenitor mass. One 
possibility is a shallower (top-heavy) IMF in the early phases 
to raise the production of massive stars, allowing progenitors 
near the black-hole mass limit (20 Mq?) to dominate enrich- 
ment of the ISM. 10 Mq SN II progenitors would neverthe- 
less still evolve within ~ 30 Myr, which is the timescale for 
the enrichment of the halo to [Fe/H] =-3.0 according to models 
which explain the rise of the r-process enrichment in [Ba/Fe] 
and [Eu/Fe] in the halo in terms of low mass (8-10 Mq) SN II 
(Mathews & Cowan [T990l Travaglio et al. 120011 Ishimaru et 
g/. .2004.1 . 
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On balance, we conclude that cuiTently there is no support 
for the proposition that only a narrow range of SN II progen- 
itor masses around 18-20 Mq contributes significantly to the 
Galactic enrichment in Mg. 

5.3.3. SN II yields 

The iron yields from SN II are notoriously difficult to pre- 
dict (Nomoto et al. I1997> : Thielemann et al. (1996 1 assume 
a decline of Fe with rising progenitor mass, while Woosley 
& Weaver J1996> assume an increase. Given this uncertainty, 
could the majority of SN II produce the same [Mg/Fe] ratio 
in their ejecta? A Mg/Fe yield ratio constant with progenitor 
mass (and thus constant over the history of the ISM) would 
require a proportionality between the size of the ejected Si- 
burned shell (the source of Fe) and the mass of the progenitor 
(which correlates with Mg yield). It seems highly unlikely that 
nature should achieve such a perfect balance in the hydrostatic 
production of Mg and the subsequent ejection of a fraction of 
the explosive production of Fe, across the ~ 10-20 or 10-30 Mq 
mass range of SN II progenitors. Moreover, the flat slope found 
for [Ca/Fe] versus metallicity would be affected by a change in 
the Fe yields, as would other element [el/Fe] ratios, and the Ca 
yield would also depend on explosion energy and progenitor 
mass. Thielemann's models do not show a variation of Ca with 
progenitor mass of the same order as that of Mg. 

Attempts to reproduce the chemical pattern of very low- 
metallicity stars by adjusting the mass cut in SN II explo- 
sion models have been unsuccessful (Chieffi & Limongi'2002'). 
Tsujimoto & Shigeyama (1998 1 traced yield versus mass for 
SN II from the observations of McWilham et al. (I1995> and 
the yields from Woosley & Weaver (1996^ and satisfactorily 
produced predictions for averaged abundances ratios within 0. 1 
dex of the observations. Their final remark is the need of revis- 
ing SN II yields of elements except C, O, and Mg. Our results 
strongly indicate that either also Mg yields are wrong (see also 
Argast et al. |2002 ), and Fe yields are wrong, or the theoretical 
framework of iGCE has reached its limit and needs a deeper 
understanding of mixing. 

Franfois et al. ( I2004> have shown that constraints to SN II 
yields can be found by matching observations with their ho- 
mogeneous GCE model. However, even if SN II yields are re- 
vised following their prescription, within the framework of ho- 
mogeneous models, the lack of scatter still remains to be ex- 
plained. If [Mg/Fe] yields truly are not constant for all progen- 
itor masses, and there is not only a very narrow mass range of 
SN II progenitors that produce significant Mg, then the lack of 
scatter points to the following unexpected conclusion: the ISM 
from which stars formed at [Fe/H] > -3.5 was on the whole 
already well-mixed, and it was enriched by a sufficient number 
of SN II to show an IMF averaged pattern at [Fe/H]< -3.0, 
in spite of the expectations to the contrary from iGCE current 
models. Furthermore, the hypothesis by Cayrel et al. f2003) of 
the existence of a plateau in the [el/Fe] ratios at very low metal- 
licity is clearly supported by our [Mg/Fe] ratio. We are left to 
puzzle over how this IMF average is achieved ove much of the 
halo on the short timescales (< 30 Myr) for the ejection of the 



full range of SN II progenitor masses (M > lOM©). This could 
imply that we are observing the primordial abundances from 
the first stars below [Fe/H] =-3.0, as Cayrel et al. suggested. 

5.3.4. Mixing and cooling timescales 

Enrichment timescales are not well constrained. Prantzos 
(.2003 ) shows that an early phase of infall and relaxing the in- 
stantaneous recycling approximation causes his halo outflow 
model to reproduce better the metallicity distribution of stars 
in the halo. Enrichment occurs much quicker in his model than 
in iGCE models, but the introduction of an early infall phase 
in his model slows the metallicity enrichment again, as a result 
of which the time needed to reach [Fe/H]=-3.0 increases from 
35 Myr to 100 Myr. A similar delay in the enrichment could 
occur in iGCE models for a strong early infall phase. However, 
a slower enrichment of the Galaxy does not affect the predicted 
scatter in [Mg/Fe] vs. [Fe/H], but rather the age-metallicity re- 
lation. 

De Avillez & Mac Low J2002t studied the mixing timescale 
in a SNe-driven ISM, and show that inhomogeneities due to 
single SN II explosions take up to 350 Myr to be erased (for 
the present Galactic SN II rate). Even when mixing scales as 
small as a few kpc are considered, efficient mixing requires of 
order 120 Myr Increasing the SN rate by a factor 10, they find 
that inhomogeneities take only some tens of Myr to disappear. 
If increasing the SN II rate reduces the time needed for effi- 
cient mixing, then for a short time (~ 30Myr) it favours the 
explosion of higher mass progenitors (because of their shorter 
lifetimes) and a higher average [Mg/Fe] ratio. This effect could 
also be achieved with a top-heavy IMF (discussed above). This 
would avoid the ISM having regions with low [Mg/Fe] values, 
but would still require a progenitor upper mass limit of 20 Mq 
to avoid the production of [Mg/Fe] > 0.4 except for at the very 
lowest metallicities (Norris (2003^)). 

We asked above how an IMF-averaged [Mg/Fe] value could 
be achieved on the very short evolutionary timescales of stars 
at [Fe/H]< -3.0. This may be the wrong question to ask if 
cooling timescales crucial to star formation, rather than mixing 
timescales, are the critical factor. iGCE models lack reliable 
treatments of the cooling of hot SN II ejecta in a metal-poor 
environment. If the cooling time exceeds the evolutionary time 
of the lowest mass SN II progenitors of 10 Mq, i.e. 30 Myr, then 
a complete IMF-average of the [Mg/Fe] ratio will be achieved, 
not on the evolutionary time of the SN II, but on the cooling 
timescale of the ISM, and therefore on the formation timescale 
of the next stellar generation. 

Oey (2003 1 investigates mixing and cooling processes in 
the ISM by means of her simple inhomogeneous model (SIM) 
of Galactic chemical evolution, which explicitly incorporates 
interstellar mixing and mass transport. As previously found by 
De Avillez & Mac Low ( I2002> . diffusion is inefficient com- 
pared to turbulent mixing. Even though "turbulent processes 
are extremely difficult to constrain", turbulent mixing is found 
to be extremely efficient in the case of a hot ionized medium 
(HIM). Thus, in a phase of HIM-dominated ISM, efficient mix- 
ing could indeed occur. Furthermore, a hot ISM has a low cool- 
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ing efficiency, thus a two-phase ISM could experience a delay 
due to the cooling time between the SN II events and the incor- 
poration of their ejecta in star forming regions. Once new stars 
could form from SN-enriched gas, they would form out of a 
well- mixed ISM. By a comparison with the observed metallic- 
ity distribution function and the lowest metallicity observed in 
stars ([Fe/H]=^.0 at the time of Oey's analysis), Oey's points 
out that an inhomogeneous phase could have been extremely 
short lived and that the ISM would have then mixed efficiently. 
Furthermore, Recchi et al. (2001 1 find that, in the case of dwarf 
galaxies, a single sturbust leads to the development of galactic 
winds and to a quick mixing. However, dwarf galaxies evolve 
in a low gravitational potential and thus dynamical processes 
(such as galactic wind) may have larger influences on the mix- 
ing than in the case of the Galaxy. 

We believe further investigations of mixing and cooling 
timescales in relation to iGCE models are advisable, since there 
is no clear theoretical expectation that the problem of the miss- 
ing inhomogeneities will be solved by adjustments of the SN II 
yields, the SN II progenitor mass range, or modification of the 
IMF. 

6. Summary 

We have presented the results of the chemical analysis of a sam- 
ple of 23 main-sequence turnoff" halo stars previously selected 
by RNB. Because of the small number of lines identified in the 
most metal-poor stars, the [Mg/Fe] abundance ratio was deter- 
mined on the basis of photometric calibrations for both effec- 
tive temperature and surface gravity. This method should avoid 
introducing star-to-star differences and preserve the homogene- 
ity of sample. 

Our results show that [Mg/Fe] is constant over the metal- 
licity range considered. There is no indication of internal scat- 
ter, the 0.06 dex standard deviation about the mean being well 
within the expected 0.08 dex random uncertainty. The maxi- 
mum point-to-point range of [Mg/Fe] is only 0.2 dex. 

The observed lack of scatter presents a challenge for iGCE 
models that predict a dispersion for [Mg/Fe] at low metallicity, 
with standard deviation of the order 0.4 dex, and a total range 
as large as 1 .0 dex. We review the possible sources of uncer- 
tainties in the iGCE models, such as the mass range of SN II 
progenitors and the Mg, Ca and Fe yields, and the IMF, but find 
no obvious solution to the problem. We identify the need for ex- 
ploring the time dependence of cooling and mixing in the early 
Galaxy, to verify whether these processes could cancel out the 
predicted scatter by ensuring that new stars form only once the 
enrichment of the ISM has been averaged over the IMF for the 
full range of SN II progenitor masses. This could occur dur- 
ing a hot early phase of the ISM, where turbulent mixing is 
extremely efficient in slowly cooling, hot, metal-deficient gas. 

It is still difficult to reconcile the large spread in [Sr/Fe] at 
[Fe/H] < -2.0 (discussed in Section 1) with the lack of a spread 
in [Mg/Fe] found in Section 4. The latter argues for star for- 
mation occurring only after an IMF average of Mg yields has 
been achieved in the ISM, whereas the lack of uniformity in 
[Sr/Fe] requires the existence of stochastic events that are not 
completely averaged out. If the source of Sr is more restricted 



in mass, or otherwise rarer than the source of Mg (which comes 
from the entire SN II range), then this apparent inconsistency 
might be resolved. The rise in [Ba/Fe] toward [Fe/H] = -2.5, 
which may be due to the origin of r-process elements in 8- 
10 Mq stars (Section 5.3.2) again argues against an IMF aver- 
age being achieved by [Fe/H] = -3.0. However, if the Mg yield 
of such low mass stars is almost insignificant, as theory sug- 
gests, then the IMF average we infer from the Mg observations 
is probably a constraint only on higher mass SN II progenitors. 
(That is to say, the IMF average claimed for Mg is an average 
over most of the IMF at masses M ~ 20 M© rather than an aver- 
age over strictly all of the SN II progenitor range.) The differ- 
ent behaviour of [Ba/Fe] and [Mg/Fe] emphasises, according 
to this interpretation, that Ba and Mg are predominantly syn- 
thesized in different progenitor mass ranges. 
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